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ABSTRACT 

Context. A segment of the early-type galaxy population hosts a prominent dust lane, often decoupled from its stellar body. Methods 
of quantifying the dust content of these systems based on optical imaging data usually yield dust masses that are an order of magnitude 
lower than dust masses derived from the observed far IR (FIR) emission. The discrepancy is often explained by invoking a diffuse 
dust component that is hard to trace in the UV or optical. 

Aims. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/sub-mm observations from the 
Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods of determining the dust content in galaxies and explore 
new ones. NGC 4370 is an edge-on, early-type galaxy with a conspicuous dust lane and regular morphology, making it suitable for 
several (semi-)analytical modelling techniques. We aim to derive the dust mass from both optical and FIR data and to investigate the 
need to invoke a putative diffuse dust component. 

Methods. We used different methods to determine the total dust mass in the dust lane. We used our exquisite optical data to create 
colour and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses were 
also derived from SED fits to FIR to sub-mm observations. Finally, inverse radiative transfer fitting was performed to investigate more 
complex dust geometries, such as an exponential dust disc and a dust ring and to treat the dust-starlight interaction in a self-consistent 
way. 

Results. We find that the empirical methods applied to the optical data yield lower limits of 3.4 x 10 5 M o , an order of magnitude 
below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses that are slightly lower, 
but fully consistent with the FIR-derived mass. We find that the effect of a nuclear stellar disc on the derivation of the total dust mass 
is minor. 

Conclusions. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a 
simple foreground screen. Moreover, by using inverse radiative transfer fitting, we are able to constrain most of the parameters that 
describe these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/sub-mm wavelengths, 
so that no diffuse dust component needs to be invoked. We furthermore caution against interpreting dust masses and optical depths 
based on optical data alone, when using overly simplistic star-dust geometries and the neglect of scattering effects. 

Key words, galaxies: individual: NGC4370 - galaxies: ISM - infrared: ISM - galaxies: fundamental: parameters - dust, extinction - 
methods: observational 


1. Introduction 

Early-type galaxies (ETGs) are not red and dead as was once the 
common understanding. Deeper observations across the electro¬ 


magnetic spectrum have revealed significant amounts of matter 
in the interstellar medium (ISM) and even star forming activity 
in these objects (see e.g. Rampazzo et al. 2005; Combes et al. 
2007; Young et al. 2011). The likely mechanism to form early- 
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type galaxies is found in major or minor mergers of spiral galax¬ 
ies. During these events, interstellar matter of the progenitors is 
partly expelled and partly consumed in waves of fast star forma¬ 
tion (see e.g. Barger et al. 1996; De Lucia et al. 2006). After a 
few Gyr, the galaxy relaxes to become a stable, non-star-forming 
object with a depleted ISM. Many ETGs, especially the more 
massive ones, are embedded in a halo of hot, X-ray-emitting gas 
(Forman et al. 1985; Sarzi et al. 2013), which is the remnant of 
the violent merger processes. Interstellar dust cannot survive in 
the harsh environment of the hot gas. Measurements of signifi¬ 
cant amounts of dust of various morphologies in over half of the 
ETGs (Goudfrooij et al. 1994c,a; van Dokkum & Franx 1995; 
Ferrari et al. 1999; Tran et al. 2001; Patil et al. 2007; Finkelman 
et al. 2008; Kaviraj et al. 2012; Kulkami et al. 2014) then lead 
to the conclusion that the dust reservoir must be replenished in 
some way. Since internal production mechanisms seem to fall 
short, a merger origin of the dust seems the most likely option 
(see e.g. Finkelman et al. 2012, for an elaborate account). 

The dust mass determination in ETGs is complicated by 
the discrepancy found between the different dust mass measure¬ 
ments. Numerous FIR observations of ETGs have been obtained 
using IRAS, ISO , and Spitzer (see e.g. Knapp et al. 1989; Leeuw 
et al. 2004; Temi et al. 2004, 2007; Xilouris et al. 2004). Temi 
et al. (2004) found that dust masses derived from ISO observa¬ 
tions (2.5 - 240 pm) are more than an order of magnitude higher 
than /RAS-derived (12 - 100 pm) dust masses. They argue that 
ISO detects much more cold dust, which could explain the dif¬ 
ference. Indeed, IRAS was only sensitive to wavelengths shorter 
than 100 pm, making it insensitive to the emission from cold 
(< 25 K) dust (Devereux & Young 1990). Because cold dust 
is actually by far the most massive dust component in galaxies, 
/RAS-derived dust masses consequently yield lower limits (Gor¬ 
don et al. 2010). 

/RAS-derived dust masses are in turn an order of magnitude 
higher than those computed using extinction models (see e.g. 
Patil et al. 2007; Finkelman et al. 2012, and references therein). 
Deriving dust masses from ultraviolet (UV), optical, and near- 
infrared (NIR) data can be difficult due to the effects of scatter¬ 
ing and the complicated star-dust distribution (see e.g. Witt et al. 
1992; Baes & Dejonghe 2001). The estimated amount of dust 
is thus highly model dependent. Moreover, it requires a consis¬ 
tent dataset spanning a sufficient wavelength range to infer the 
wavelength-dependent extinction by dust. Up to now, the gener¬ 
ally applied dust geometry was that of a dust screen in between 
the galaxy and the observer (see e.g. Goudfrooij et al. 1994a; 
Ferrarese et al. 2006; Patil et al. 2007; Finkelman et al. 2008, 
2010). This naturally yields lower limits of the dust mass when 
derived from extinction. The dust mass discrepancy with far-IR 
(FIR) measurements may be explained by invoking a second dust 
component, which follows the stellar body and is diffuse (hence 
not detectable in the optical; Goudfrooij & de Jong 1995; Wise 
& Silva 1996). 

The diffuse dust component in ETGs is not to be confused 
with the component of the same name found in spiral galaxies. 
The latter term was invoked to differentiate the dust that resides 
in clumps or filaments in the disks of spiral galaxies. The dif¬ 
fuse dust in ETGs does not refer to this low-scale height compo¬ 
nent, but to the dust components following the stellar distribu¬ 
tion. Furthermore, in ETGs, it is heated solely by the old stel¬ 
lar populations, while in the discs of spiral galaxies, young stars 
contribute significantly to the heating of the diffuse dust (see e.g. 
De Looze et al. 2014; Bendo et al. 2015, and references therein). 

To explain the dust mass discrepancy found in ETGs, a dif¬ 
fuse dust component should be about 10-100 times as massive 


as the dust lane (based on IRAS and ISO mass estimates, re¬ 
spectively). Such a vast reservoir of cold dust should then be 
observable with the Herschel Space Observatory (Pilbratt et al. 
2010). However, in most ETGs, the FIR and sub-mm emission 
is much more compact (see e.g. Smith et al. 2012b). As a result, 
dust emission is not spatially resolved in the vertical direction 
in Herschel observations. Therefore, a direct detection of a dif¬ 
fuse dust component following the stellar distribution has not yet 
been reported in ETGs. 

We can now revisit these problems using high-quality data 
covering a broad range of wavelengths and advanced models. 
In the UV, optical and NIR, large surveys such as SDSS (York 
et al. 2000), 2MASS (Skrutskie et al. 2006), GALEX (Martin 
et al. 2005), and WISE (Wright et al. 2010) have made multi¬ 
wavelength datasets available for many objects. In the FIR/sub- 
mm regime, Herschel proved a tremendous leap forward in terms 
of sensitivity. Consequently, the dust content in ETGs has been 
studied intensely over the past few years. Dedicated investi¬ 
gations of individual targets (e.g. M87 Baes et al. 2010, M86 
Gomez et al. 2010, and NGC 4125 Wilson et al. 2013) were 
supplemented with statistical studies. For example, from the 
Herschel Reference Survey (HRS, Boselli et al. 2010): Smith 
et al. (2012b); Cortese et al. (2012b), from the Herschel Virgo 
Cluster Survey (HeViCS, Davies et al. 2010): Clemens et al. 
(2010); di Serego Alighieri et al. (2013), and from the Herschel 
Astrophysical Terahertz Large Area Survey (H-ATLAS, Eales 
et al. 2010): Agius et al. (2013); Rowlands et al. (2012). As the 
SPIRE (Griffin et al. 2010) and PACS (Poglitsch et al. 2010) in¬ 
struments on-board the Herschel spacecraft cover wavelengths 
from 70-500 pm, the peak of cold dust emission could be suffi¬ 
ciently sampled, and more reliable dust masses could be derived 
for all of these objects. 

On the modelling side, significant progress has been made in 
the field of dust radiative transfer (for a review, see Steinacker 
et al. 2013). Multiple efforts have already succeeded in fitting 
radiative transfer models to individual galaxies (Xilouris et al. 
1997, 1999; Popescu et al. 2000; Bianchi 2007; De Looze et al. 
2012b,a; De Geyter et al. 2013, 2014). Up to now, this has 
mainly been applied to edge-on spiral galaxies, which inherently 
have morphologies that are more complex than ETGs. A similar 
way of modelling can thus be applied to ETGs with a prominent 
dust lane. 

In this paper, we re-evaluate the dust mass discrepancy in 
ETGs using new high-quality optical and FIR/submm data, and 
advanced models. NGC 4370 is an edge-on lenticular in the 
Virgo cluster at a distance of 23 Mpc (taken from the GoldMine 
archive; Gavazzi et al. 2003). It has been classified as a boxy SO 
galaxy in several independent studies (Binggeli et al. 1985; Bail- 
lard et al. 2011; Buta et al. 2015). However, the edge-on view 
makes it difficult to classify because stellar discs or spiral arms 
may be hidden by the dust lane. A difference of more than a fac¬ 
tor 3 between the dust mass obtained from the modelling of opti¬ 
cal data and from IRAS fluxes was found for this galaxy (Finkel¬ 
man et al. 2008). Since IRAS dust masses are usually underesti¬ 
mations, this discrepancy can only widen. Thanks to its smooth 
stellar profile and conspicuous, regular dust lane, the galaxy is 
suited to several (semi-)analytical modelling techniques. Excel¬ 
lent optical data are now available for this object, complemented 
by deep FIR/sub-mm observations from Herschel. We therefore 
choose NGC 4370 for this pilot study. The paper is organised 
in the following way. The data is described in Sect. 2, and our 
modelling efforts are outlined in Sect. 3. We discuss the results 
in Sect. 4 and present our conclusions in Sect. 5. 
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2. Data 

2.1. Optical/NIR data 

Optical imaging for NGC 4370 was obtained from the Next 
Generation Virgo Cluster Survey (NGVS, Ferrarese et al. 2012). 
This legacy survey was executed with the MegaCam wide-held 
optical imager at the Canadian-France-Hawaii-telescope over 
the 2009A-2013A cycles. All fields were observed in the ugiz 
filters, with some also covered in the r band. The survey con¬ 
sists of 117 separate MegaCam fields of - 1.0 deg 2 with a pixel 
scale of 0.187 arcsec. The fields overlap slightly, covering a total 
area of more than 104 deg 2 on the sky. The fields are centred on 
M87 and M49, and they map the Virgo cluster out to one virial 
radius. The NGVS data is photometrically calibrated using stars 
from the SDSS resulting in rms errors of 0.01 - 0.02 mag (Gwyn 
2008). Through dedicated data processing, a point-source depth 
of ~ 25.9 mag (10cr) is achieved in the g band. For extended 
sources, the surface brightness limit is g = 29 mag arcsec -2 for 
a 2 cr detection above the sky. 

NGC 4370 was imaged in the ugi and z bands. This dataset 
exposes the galaxy in great detail (seeing < 1") and allows for an 
accurate study of the extinction effects (see Fig. 1). The fully re¬ 
duced and calibrated images were first brought to the same point 
spread function (PSF). The full width half maximum (FWHM) 
of the PSF was determined for each band as an average of sev¬ 
eral brighter stars in the image. The resolution (FWHM) was 
measured to be 0.81, 0.87, 0.53, and 0.70 in the u , g , i, and z 
bands, respectively, and therefore all images were convolved to 
a common 0.87 arcsec seeing using a Gaussian convolution ker¬ 
nel. The frames were regridded to the pixel grid of the g band as 
well in order to work with a set of pixels that correspond to the 
exact same physical regions. After that, the images were rotated 
over 6° in such a way that the dust lane fell along the x-axis of 
the frame. Both steps required some interpolations that reduced 
our working resolution to 0.98". This corresponds to 97 parsecs 
at the distance adopted for NGC 4370 (23 Mpc). 

In a final step, bright foreground stars were masked out in 
each band. The stars were identified based on their shape and 
offset to the galaxy’s brightness profile and verified with the NO¬ 
MAD catalogue (Zacharias et al. 2004). Background galaxies in 
the field of view were identified using NED 1 . All reported galax¬ 
ies in the field lie near the edge of the image, and the brightest 
one has a g band magnitude of only 20.3. We therefore did not 
mask any background sources. 

The dust in NGC 4370 greatly affects the appearance of the 
optical images, which we use to our advantage. However, it is 
also useful to gain more insight into the underlying stellar mor¬ 
phology. We therefore include the IRAC 3.6 yum image from the 
Spitzer Survey of Stellar Structure of Galaxies (S 4 G, Sheth et al. 
2010) which is shown in the top right-hand panel of Fig. 1. At 
these wavelengths, we expect minimal dust contamination and 
still benefit from a decent spatial resolution (FWHM ~1.7 arc¬ 
sec). 

2.2. FIR/sub-mm data 

The Herschel Virgo Cluster Survey (HeViCS; Davies et al. 2010) 
observed 84 deg 2 of the Virgo Cluster with both PACS (100, 160 
/rm) and SPIRE (250, 350, 500 yum). The observations were done 
in parallel fast-scan mode (60 arcsec s -1 ). The area is subdivided 
into four areas, corresponding to the region around M87 (V2), 
and the clouds to the north-west (VI), west (V3), and south (V4) 

1 http://ned.ipac.caltech.edu 
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Fig. 1: Overview of the dataset. Top left: RGB colour map of NGC 
4370 based on the giz images of the NGVS. Top right: IRAC 3.6 yt/m 
image from the S 4 G. Bottom: PACS 100 yum (left) and 160 yum (right) 
images. The beam size is indicated as a green circle in the lower left 
comer of each image. 


of V2. Each area was covered by a 4 x 4 sq. degree tile and 
observed four times in two orthogonal scan directions to reduce 
the 1 // noise. The PACS data were reduced using a combina¬ 
tion of HIPE (vl2.0, Ott 2010) and Scanamorphos (version 24, 
Roussel 2013). The latter, in particular, corrects for a Hatfield 
distortion that causes an offset of up to -7% in the fluxes at 160 
/im. In a subsequent observing run, 44 galaxies from the bright 
galaxy sample of HeViCS were observed with PACS in the 70 
and 160 pm bands in map-making mode. NGC 4370 was one 
of them, and we use the 70 pm data from this run. These obser¬ 
vations were reduced in a similar way to the parallel scan maps. 
The SPIRE observations were reduced using HIPE (vl2.0) and 
the custom routine BriGAdE (Smith, in prep.) For more details 
on the data reduction, we refer the reader to Auld et al. (2013). 
The final maps have a beam size (FWHM) of 5.7, 9.4, 13.4, 
18.2, 25.4, and 36.0 arcsecs for 70, 100, 160, 250, 350, and 500 
yum, respectively, corresponding to a physical scale of 0.64, 1.05, 
1.49, 2.03, 2.83, and 2.01 kpc at a distance of 23 Mpc for NGC 
4370. 

Because NGC 4370 lies on the intersection between two tiles 
(V3 and V4), it was observed -partially- by scans targeting the 
two different fields. It turns out that the galaxy was observed 
by a total of 13 scans, which were thus used to construct the 
maps at 100 and 160 //m. This is allowed by the option cutout 
in Scanamorphos, which forces the use of only those parts of 
the scanlegs that contain emission from the galaxy. This will en¬ 
hance the signal-to-noise ratio (S/N), positively affecting the flux 
uncertainty determination, so we chose to redo the flux density 
measurement in the PACS 100 pm and 160 yum bands. We find 
new values of Fiooyum = 3.03 ± 0.18 Jy and F\ 6 0/im = 3.65 ± 0.24 
Jy. This is about 3% lower for Fioo/im and 6% for F i 60 /im 
compared to the measurements of Cortese et al. (2014), who 
used the standard data products. These are, however, not cor¬ 
rected for the Hatfield distortion in the PACS bands, while our 
measurements do include this correction. Figure 1 shows the 
newly reduced PACS maps. We also measure the integrated flux 
of the PACS 70 pm image and find F 10yU m = 1.348 ± 0.073 
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Jy, which is fully consistent with the MIPS 70 p m flux mea¬ 
sured by Bendo et al. (2012). Since the SPIRE flux measure¬ 
ments from Ciesla et al. (2012) are already of good quality with 
just eight cross-scans, we take these for our further analysis. 
They are F 25011 m = 1.95 ± 0.05, T 35 o^m = 0.77 ± 0.04, and 
^ 500 yum = 0-25 + 0.02 Jy. Table 1 summarises the integrated 
fluxes for NGC 4370 used in this work. 


3. Methods and results 

We explore several methods of investigating the extinction prop¬ 
erties and determine the dust mass of NGC 4370. Each model 
treats the dust geometry, hence the corresponding attenuation ef¬ 
fects, in a different way. It is important to note that different 
dust mass estimates cannot be compared unless a common dust 
model is used. To achieve this consistency, we adopt the Draine 
& Li (2007) dust model throughout the paper. The dust model is 
derived from observations in our own Galaxy and consists of car¬ 
bonaceous grains, amorphous silicates, and polycyclic aromatic 
hydrocarbons (PAH). 

3.1. Dust masses from FIR/sub-mm emission 

The bulk of the dust in galaxies resides in the diffuse ISM at low 
equilibrium temperatures of 10 - 30 K. This component is di¬ 
rectly measurable through its emission in the FIR/sub-mm. The 
emission is mostly optically thin and so guarantees that all the 
dust will be visible, unlike in the case of the UV/optical/NIR 
where some dust may be hard to trace. Although the SPIRE ob¬ 
servations were confusion limited, it may still be that we miss 
faint emission at larger radii. From the stacking of HRS spi¬ 
ral galaxies, Smith et al (in prep.) find that ~2.4% of the dust 
resides outside the D 25 radius. Assuming this is even less for 
ETGs, the missing dust will not contribute significantly to the 
sub-mm emission. 

In the past, the dust mass and temperature of dusty ETGs 
were derived using a modified black body function (Hildebrand 
1983) and IRAS fluxes (e.g. Young et al. 1989; Goudfrooij & de 
Jong 1995; Patil et al. 2007; Finkelman et al. 2008). For NGC 
4370, in particular, a dust mass of 4.54^ ^10 5 M© was found, 
scaled to the distance we assume in this work (23 Mpc). The dust 
temperature from the IRAS colours was found to be T& = 43.1 ± 
1.3 K (Finkelman et al. 2008). It must be noted, however, that 
IRAS was unable to detect dust with temperatures colder than 25 
K. The values thus only apply to the warm dust component. 

Since the broad availability of Herschel observations, it has 
become common practice to estimate the dust mass from a mod¬ 
ified black body fit to the FIR and sub-mm fluxes. Following 
the recommendations of Bianchi (2013), we only fit a modified 
black body model with a fixed J3 = 2. This allows us to deter¬ 
mine the dust mass using the Draine & Li (2007) dust model, 
with an emissivity of K^o/im = 0.192 m 2 kg -1 (see Fig. 2). We 
find a temperature T d = 21.8 ± 0.3 K and a total dust mass 
^modBB _ 6J7+0.38 x i () 6 M 0 . This is in line with di Serego 
Alighieri et al. (2013), who used similar data and methodology. 
It is roughly an order of magnitude higher than the TRAS-derived 
dust mass and states again the importance of sampling the sub- 
mm emission of the cold dust component. 

We also tested a more advanced spectral energy distribution 
(SED) fitting model: MAGPHYS (da Cunha et al. 2008). This 
code allows simultaneous treatment of both stellar and dust com¬ 
ponents, preserving the energy balance between them. It has 
been used on a wide variety of galaxies (e.g. da Cunha et al. 


Table 1: Observed fluxes used for the SED fits for NGC 4370. 


Band 

Flux (mJy) 

Source 

FUV 

0.0463 ± 0.0063 

a 

NUV 

0.2166 ± 0.0059 

a 

u 

6.1 ± 1.1 

b 

g 

22.9+1.8 

b 

i 

62.3 ± 2.5 

b 

z 

80.9 ± 3.2 

b 

J 

107 + 21 

c 

H 

138 + 21 

c 

K 

121 + 24 

c 

WISE 1 

73.8 + 1.9 

d 

WISE 2 

39.8 + 1.3 

d 

WISE 3 

64.8 + 3.1 

d 

WISE 4 

56.4 + 4.7 

d 

MIPS 24 pm 

48.3 + 2.1 

e 

MIPS 70 pm 

1284+ 132 

e 

PACS 70 pm 

1348 + 73 

f 

PACS 100 pm 

3027 + 172 

f 

PACS 160 pm 

3647 + 231 

f 

MIPS 160 pm 

2880 + 350 

e 

SPIRE 250 pm 

1947 + 53 

f 

SPIRE 350 pm 

772 + 38 

f 

SPIRE 500 pm 

251 + 19 

f 


References, (a) Cortese et al. (2012a), (b) NGVS, this work, (c) Gold- 
Mine (Gavazzi et al. 2003), (d) Agius et al. (2015), (e) Bendo et al. 
(2012), (f) HeViCS, this work. 


2010; Smith et al. 2012a; Clemens et al. 2013; Lanz et al. 2013; 
Viaene et al. 2014; Baes et al. 2014). We sampled the panchro¬ 
matic SED of NGC 4370 using GALEX fluxes from Cortese et al. 
(2012a), the ugiz fluxes from the NGVS, JHK magnitudes from 
the GoldMine archive (Gavazzi et al. 2003), WISE fluxes from 
Agius et al. (2015), and MIPS fluxes from Bendo et al. (2012). 
The Herschel data described above complete our dataset. Table 1 
lists the fluxes that were used. The best-fitting SED is shown in 
Fig. 2, as is the intrinsic (stellar) SED without dust attenuation. 
The model is able to match the observed data remarkably well. 
If we use the same Draine & Li (2007) dust model, we find a to¬ 
tal dust mass of 5.7^^ x 10 6 M© and a temperature of 22A + _^\ 
K for the cold dust. It is reassuring that we retrieve -within the 
error bars- a similar cold dust temperature and total dust mass to 
the modified black body model. 

We note that even these dust mass determinations are still 
only a best guess of the true dust mass. There are known caveats 
in the modified black body fitting method. First, a single tem¬ 
perature is assumed for the cold, diffuse dust. In reality, this is 
likely a continuous range of temperatures, heated by a combina¬ 
tion of emission sources. The multiple thermal components of 
dust appear to blend in such a way that it can be fit by a single 
modified black body (Bendo et al. 2015). In this scenario, the 
resulting dust temperatures from a single modified black body 
may be higher than the actual dust temperatures, and the result¬ 
ing dust masses may be biased to lower values. Second, it is 
assumed that the dust in NGC 4370 has the same mixture as in 
our own Galaxy and that this mixture is the same everywhere in 
the galaxy. This should not, however, affect relative dust mass 
differences such as those investigated in this paper. Last, recent 
results from Plank Collaboration et al. (2014) indicate that the 
Draine & Li (2007) dust model overestimates the dust mass by a 
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Fig. 2: MAGPHYS fit of the panchromatic SED of NGC 4370. The 
black points are the observed fluxes, the blue line is the unattenuated 
(stellar) SED, and the red line the attenuated, best-fitting template SED. 
The residuals are plotted below. The dashed green line is the best-fitting 
modified black body function to the Herschel fluxes. 


factor 2-4. Throughout this work we nevertheless consider SED 
fitting to be the most objective and reliable way to estimate at 
least the order of the true dust mass in NGC 4370. 

3.2. Dust masses from colour maps 

Given the homogeneous dataset, colour maps are the simplest 
way to trace the extinction by dust. We can write the observed 
colour between two bands X and Y as 


(X-Y) = (X 0 -Y 0 ) + (A x -A y ) (1) 

where Xq - Yq is the intrinsic colour and A x -A Y the difference in 
attenuation between the two wavelength bands, also called selec¬ 
tive attenuation. The intrinsic colour Xq - Yq is a manifestation 
of the wavelength-dependent luminosity of the stellar popula¬ 
tions. The intrinsic colour was determined from the regions in 
the galaxy that do not show any obscuration by dust. Because we 
have four bands, six colour maps can be constructed: u-g,u-i , 
u-z, g-i, g~z, and i-z. Figure 3 shows radial profiles along the 
major and minor axes for each of the colours. The same profile 
is retrieved in all colours: along the major axis, reddening starts 
increasing around ~30 arcsec from the centre when dust redden¬ 
ing starts to be effective. It peaks in the centre where most dust 
is expected along the line of sight. Along the minor axis, a small 
gradient is present outside the dust lane. The reddening increases 
sharply when entering the dust lane at ~4 arcsecs and peaks in 
the centre. Not surprisingly, the highest reddening is found for 
the u - z colour and the lowest reddening for i - z . 

Interestingly, at the centre of the dust lane, a minor dip is 
observed in the colour profiles along the minor axis. This effect 
is most evident in the u-g profile (dark green), but also visible in 
most of the other profiles. Figure 4 shows the maps of selective 
attenuation, i.e. colour maps minus intrinsic colour. Here, the 
central part of the dust lane also shows this dip in reddening as 
a dark line running along the major axis. This colour inversion 
effect has been observed by Bianchi (1999) in the B-I colour map 
of NGC 891, a massive edge-on spiral galaxy. It is suggested that 
saturation of the dust extinction occurs in these regions. At the 
same time, the number of stars between the dust and observer is 
expected to be largest here. As a result, the colour of the stars 
that lie in front of the dust becomes dominant over the reddening 
by dust. This makes the central part appear bluer than the dusty 
regions right above or below. 


We must note that there are known colour gradients in ETGs, 
which possibly affect the measurement of the intrinsic colour. 
Petty et al. (2013) showed that the inner parts of ETGs are gen¬ 
erally redder than the outskirts, and several spectroscopic stud¬ 
ies have revealed age and metallicity gradients (e.g. Carollo et al. 
1993; Kobayashi & Arimoto 1999; Baes et al. 2007; Koleva et al. 
2011). Asa result, when measuring the intrinsic colours of NGC 
4370 outside the dust lane, it is possible that we pick up bluer 
colours than those behind the dust lane. This would result in 
slightly lower A x - Ay values, hence lower dust masses. We 
confirm colour gradients along the minor axis of NGC 4370 (see 
Fig. 3), however they are small with respect to the reddening in 
the dust lane. We therefore assume that the above effect will only 
cause a limited decrease in dust mass. The resulting maps of se¬ 
lective attenuation are shown in the left-hand column of Fig. 4. 

On the other hand, since NGC 4370 is classified as an SO type 
galaxy, a stellar disc may be present, concealed by the dust lane. 
The IRAC 3.6 pm image (see Fig. 1) indeed shows a disc-like 
feature in the centre of the galaxy. If this obscured disc contains 
some ongoing star formation, it will be bluer than the intrinsic 
colours we measure outside the dust lane. Consequently, this 
can lead to underestimating the dust mass. It is nearly impos¬ 
sible to determine the intrinsic colour of a potential stellar disc 
from the edge-on perspective, which immediately illustrates the 
limitations of this simple approach. Fortunately, the stellar disc 
is relatively small (< 25" in diametre or < O. 25 D 25 ) and lim¬ 
ited to the central regions of NGC 4370, so this effect may be 
relatively small. Nevertheless, this implies that resulting dust 
masses should be seen as lower limits to the true dust mass in 
NGC 4370. 

The degree of the reddening is limited (1-4 magnitudes) be¬ 
cause we only consider optical observations. The reason not to 
expand this set to UV or NIR observations is dual: first is the su¬ 
perb resolution in the optical images. The FWHM of the IRAC 
maps is twice that of the g band, the largest optical beam size. 
For GAFEX observations, the difference is more than a factor 
of 5. Second are the intrinsic colour gradients, which are much 
stronger when colour maps are created over these wide wave¬ 
length ranges. An accurate subtraction of the intrinsic colour is 
crucial for producing a reliable selective attenuation map. We 
therefore chose to proceed with the limited but fully consistent 
and high-quality optical dataset. 

To convert (selective) attenuation maps to dust maps, we 
need to make assumptions about the star-dust geometry. Dust 
masses can be derived from the dust mass surface density in each 
pixel Sj. The total extinction coefficient k a of the dust at a cer¬ 
tain wavelength, relates the dust mass to the optical depth t a : 


t a = K A Y, d . ( 2 ) 

To determine the optical depth from the attenuation maps, the 
geometry of the dust distribution needs to be taken into account. 
We consider two different dust geometries. The first one is a sim¬ 
ple foreground dust screen between the stars and the observer: 


A a = -2.5 lo g(e Ta ) = 1.08 6t a . 


(3) 


This geometry requires a minimal amount of dust to obtain the 
observed attenuation. The second model describes a layer of dust 
embedded in a layer of stars and is usually dubbed “the sandwich 
model” (Disney et al. 1989; Boselli et al. 2003; Cortese et al. 
2008): 


A a - -2.5 log 




(4) 
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Fig. 3: Colour profiles along the major and minor axes, derived from 
each colour map of NGC 4370. 

The dimensions of the dust disc, relative to the stellar disc, are 
expressed by the parameter We limit ourselves to the two 
extreme cases of this geometry, —> 1 and —> 0. The former 
limit corresponds to a dust layer of equal dimensions to the stel¬ 
lar layer. This boils down to a uniform mixing of dust and stars: 

A a = -2.5 log I———I. (5) 

This formula can in principle be inverted analytically by means 
of the Lambert W function to recover the optical depth from an 
observed attenuation (see Boquien et al. 2013). 

On the other hand, when —> 0, an infinitely thin slab of 
dust is embedded in the stellar layer along the line of sight: 

/1 + e~ TA \ 

A a = -2.5 log I—-—I. ( 6 ) 

Colour maps yield maps of selective attenuation A Al - A Al . 
This can be written as a combination of optical depth functions 
using equations (3) - ( 6 ): A M - A Al = f(r M ) - /(r^). Using the 
Draine & Li (2007) dust model, the optical depth at one wave¬ 
length can be converted into the equivalent optical depth at an¬ 
other wavelength: 

_ TC ^ rj ^ 

KAx KA 2 


Offset from centre (arcsec) 
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Fig. 4: Left: Selective attenuation ( A Al - A Al ) for all optical colours 
of NGC 4370. Right: corresponding dust mass maps in the foreground 
screen geometry at the critical 1 arcsec resolution. The white dots indi¬ 
cate masked foreground stars. 


These equations were then numerically solved to yield a value 
for a specific optical depth, which in turn yields a dust mass 
through Eq. 2. 

Uncertainties on the attenuation and attenuation differences 
were derived following Goudfrooij et al. (1994b). Attenuation 
values were first derived on a pixel-by-pixel basis at the native 
pixel scale of 0.19". This is roughly five times smaller than the 
PSF of the convolved images (0.98"). We therefore rebinned 
each attenuation image using a 5 x 5 pixel box and assigned to 
each (rebinned) pixel a value and uncertainty equal to the mean 
attenuation (or selective attenuation) and standard deviation cal¬ 
culated within each 5x5 pixel box. These uncertainties are prop¬ 
agated further through the standard bootstrapping method; we 
created 1000 randomly perturbed attenuation (or attenuation dif¬ 
ference) maps based on the errors of the real attenuation maps. 
For each map, the resulting dust mass is then computed. The 
uncertainty on the dust mass is then taken to be three times the 
standard deviation of the random set of dust masses. The final 
dust masses and their errors are listed in Table 2. 

Figure 4 shows the dust mass maps derived from the corre¬ 
sponding colour maps in the case of the foreground screen ge¬ 
ometry. We only sum pixels inside the dust lane, because this is 
the only region where dust is directly visible. In the foreground 
dust screen geometry, the same morphology is retrieved for all 
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dust mass maps. Although there are some local variations in the 
dust mass per pixel, the total dust masses lie close to each other 
and average out to 5.60 ;^^ x 10 5 M Q , which should be inter¬ 
preted as a lower limit. This is a value that is comparable to the 
/RAS-derived dust mass, but still an order of magnitude below 
the dust mass estimates using Herschel data. 

In the case of the sandwich model, retrieving a real solution 
for the optical depth was not possible for many pixels. For the 
u-z colour, for example, the maximal attainable selective atten¬ 
uation based on eq. (5) is 1.2597. For eq. (6), this is only 0.3276. 
In contrast, the observed selective attenuation for the u -z colour 
reached 2.49. Consequently, the pixels with the largest selective 
attenuation are incompatible both with an embedded thin slab 
of dust (£i —> 0) and with a uniform mixing of dust and stars 
(£i —> 1). Their mixed nature means that changing the proper¬ 
ties of the dust distribution also changes the location of the stars 
relative to the dust. Because there will always be stars in front of 
and behind the dust, there are inherent limitations to the attenu¬ 
ation values along a line of sight. It appears that the sandwich 
model - in the limits we explore here - cannot achieve attenu¬ 
ation values that are high enough to explain the observed ones. 
We therefore chose not to pursue these geometries any further 
for the colour maps. 

3.3. Dust masses from multi-Gaussian expansion models 

As a third method, we compare the observations with a dust free 
model of the galaxy in each waveband to construct monochro¬ 
matic attenuation maps. A dust-free model is constructed 
through multi-Gaussian expansion (MGE) fitting. This tech¬ 
nique was first developed by Emsellem et al. (1994) and consists 
of the summing of multiple 2-D Gaussian functions to mimic the 
brightness profile of elliptical galaxies. Cappellari (2002) con¬ 
structed a dedicated set of IDL routines to efficiently fit a set of 
Gaussians to an observed image. 

It is crucial that the modelling is based on the dust-free part 
of the galaxy so that the dust lane is masked. We fixed the co¬ 
ordinates of the centre of the galaxy as determined from the z 
band image (which is assumed to have the least dust contami¬ 
nation) using the MGE-IDL routine find_galaxy. NGC 4370 
has a boxy stellar distribution, not easily modelled by an ellip¬ 
tical profile. The MGE fitting software fortunately allows each 
Gaussian to have a different position angle, relative to the x-axis 
of the frame. A boxy brightness profile can be constructed by 
summing multiple 2-D Gaussian distributions with a slightly dif¬ 
ferent angle. 

MGE models were fitted to each of the four bands of our 
data set. A roughly equal effective area was used to constrain 
the free parameters. This corresponded to a S/N threshold of 
three in the u and g bands, and five in the i and z bands. The 
number of initial Gaussians is fixed for each MGE run. Allowing 
more Gaussians increases the number of free parameters, so that 
more complex brightness profiles can be constructed. On the 
other hand, if too many free parameters are available, peculiar 
and unphysical profiles may be found. To find a balance, we 
evaluated the non-reduced^ 2 and shape of the brightness profile 
as a function of the number of free parameters until the best fit 
was obtained. 

The resulting set of Gaussians allows constructing a 2-D, 
dust-free model of the galaxy. We can now compute the attenu¬ 
ation A a for each of the bands by taking 

A, = -2.5 log (8) 

\G,model / 


where G,model and G,image are the model and reference frames, 
respectively. The attenuation values can then be converted to 
optical depths using equations (3), (5), and (6). Equation (2) 
then allows the computation of the dust mass maps. 

Figure 5 shows the results of the MGE modelling for the ugiz 
bands. In the first column, the observations are overlaid on the 
contours of the dust-free MGE model. In all bands, the contours 
closely follow the stellar emission in the observed images. Gen¬ 
erally, the model is elliptical in the centre and becomes boxier 
in the outer part of the galaxy. Owing to the obscuration of the 
dust lane, some small differences in the central contour shapes 
are visible across the wavebands. The second column shows the 
residuals of the MGE model and the observations. Obviously, 
the dust lane sticks out here because no dust was included in 
the MGE model. Apart from that, the residual levels are usually 
below 12.5% and only increase in the outer, low S/N regions. 
Using this formalism, it is thus possible to reproduce the stellar 
emission of NGC 4370, and we can now construct attenuation 
maps. These are shown in Fig. 8 and discussed in Sect. 3.4. 

In the last two columns of Fig. 5, the resulting dust mass 
maps are displayed for the screen, along with uniform-mix ge¬ 
ometries. In the case of the embedded thin slab (eq. (6)), again 
no real solutions were found for many of the pixels. This partic¬ 
ular geometry can only generate moderate attenuation values up 
to 2.5 log 2 = 0.753, which are lower than most pixels in the dust 
lane. We therefore do not pursue this geometry any further. For 
the other two geometries, the dust mass maps have similar mor¬ 
phology and pixel values. A slight decrease in pixel dust-mass 
values is observed in the z band. This is not entirely surprising 
since the z band is the longest in wavelength, hence least sensi¬ 
tive to attenuation by dust. Comparing both the dust screen (3rd 
column) and uniform dust mix (4th column) geometries, the lat¬ 
ter clearly has higher dust masses in each pixel. Owing to the 
mixing of dust and stars, more dust is needed to obtain the same 
level of attenuation as for the foreground screen of dust. The 
obtained dust masses are consistent with this view. The total 
masses are listed in Table 2 and were obtained by summing the 
pixel values in the dust lane. As already visible from the maps, 
the dust masses in the u , g , and i bands lie close to each other, 
within each geometry. The z band total dust masses are signifi¬ 
cantly lower. 

On average, the foreground screen yields a total dust mass of 
3.35 ^qX10 5 M o and the uniform mix geometry 7.9^'^xl0 5 M o . 
These values are comparable to the average total dust mass from 
the colour maps and consistent with the TRAS-based dust mass 
for this galaxy. In the end, with this high-quality data, we are 
already able to alleviate the previously found dust mass discrep¬ 
ancy between optical and /RAS-based dust masses. Of course, 
IRAS masses are still insensitive to cold dust. Indeed, as for the 
colour maps, the masses are still about an order of magnitude 
lower than the ones obtained using the Herschel observations, 
which do observe emission from cold dust. 

3.4. Dust masses from radiative transfer modelling 

The most advanced way of treating the extinction of starlight 
through a dusty medium is through radiative transfer (RT) sim¬ 
ulations. This technique allows more realistic dust distributions 
(e.g. embedded discs or rings) and takes the 3-D geometry of 
the problem into account. The main advantage over the previ¬ 
ous techniques is the self-consistent treatment of dust extinction 
along the line of sight. This provides a more realistic determi¬ 
nation of the amount of dust needed to produce the observed 
dust lane. Ideally, we would like to solve the inverse RT prob- 
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Fig. 5: First column: Observed images of NGC 4370, overlaid with contours of the corresponding dust-free MGE model. The contour levels are 
the same for all bands. 2nd column: Residuals of the MGE models and the observed images. 3rd column: corresponding dust mass maps in the 
foreground screen geometry and at the critical 1 arcsec resolution. 4th column: corresponding dust maps in the uniform star-dust mix geometry. 
The white dots indicate masked foreground stars. 


lem, i.e. find the best fitting parameters for a RT-model to match 
the observations. FitSKIRT (De Geyter et al. 2013, 2014), an 
optimization code built around the RT code SKIRT (Baes et al. 
2003,2011; Camps & Baes 2015), is able to probe this parameter 
space efficiently and find the best-fitting model to a set of obser¬ 
vations. FitSKIRT’s optimization technique is based on genetic 
algorithms, which start with an initial population of 150 differ¬ 
ent parameter sets. For each set, a RT simulation is run and the 
simulated images are compared to the observed images. Then, 
a second generation of 150 parameter sets is constructed from 
the best-fitting sets. Again, a RT simulation is run for each pa¬ 
rameter set. We chose to limit the number of generations to 100, 
although convergence is usually reached much faster. The code 
was originally designed and tested to fit the complex geometries 
of edge-on spiral galaxies. These objects are successfully mod¬ 
elled with up to 19 free parameters. Several of those describe 
the complex bulge-disc geometry of the stars. In our model, 
however, we keep the stellar geometry fixed, which drastically 
reduces the number of free parameters to fit, thereby increasing 
our confidence in the results of the fits. 


Imaging in blue bands is ideal for investigating the extinction 
properties of dust, but it hides information on the underlying stel¬ 
lar distribution. The red bands, on the other hand, provide more 
accurate information on the stellar distribution, but it makes it 
more difficult to study the dust-extinction effects. A combined 
fit to all bands can alleviate these degeneracies, as shown by De 
Geyter et al. (2014). This procedure is called oligochromatic RT 
fitting. An additional benefit is the more stable convergence to¬ 
wards the optimal solution, particularly for low S/N observations 
(HauBler et al. 2013; Vika et al. 2013). We explore RT models 
with two different geometries for the dust distribution in NGC 
4370: an exponential disc geometry and a ring geometry. Their 
best-fitting parameters are determined using FitSKIRT. We si¬ 
multaneously fit each RT model to the four images ( u , g , /, and z 
bands). 

The stellar distribution for each model is the same and de¬ 
rived from the z-band image. The longest waveband was cho¬ 
sen to limit contamination by the dust. We construct a dust-free 
MGE model similar to the one obtained in section 3.3. In this 
case, however, we force each Gaussian to have the same posi- 
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Table 2: Dust masses for NGC 4370 as derived from the various meth¬ 
ods, in units of 10 5 M Q . The embedded slab geometry did not yield real 
solutions for either the colour or the MGE method and are therefore not 
included in the table. 


Colour maps 

Screen 

Mix 

u-g 

4.48 + 0.15 

no solution 

u - i 

5.37 + 0.15 

no solution 

u- z 

5.48 + 0.14 

no solution 

g-i 

6.03 ±0.17 

no solution 

g-z 

6.06 + 0.17 

no solution 

i-z 

6.20 ± 0.27 

no solution 

Average 

5.60+°f 

no solution 

MGE fits 

Screen 

Mix 

u 

3.25 + 0.12 

8.20 + 0.31 

g 

3.83 + 0.15 

9.36 + 0.37 

i 

3.58 + 0.24 

8.14 + 0.55 

z 

2.72 + 0.27 

5.93 + 0.49 

Average 

7} Q c+0.64 
-0.90 

7 0+1.9 
' -2.5 

IRAS a 

ModBB 

MAGPHYS 

FIR/sub-mm 

4 C 4 + 0.33 
-0.31 

61 7 +3,8 

Ui -3.6 

57+f 


FitSKIRT 

Exponential disc 

50 + 24 

Ring 

33 + 27 

IRAC Ring 

70+18 


References, (a) Finkelman et al. (2008) 

tion angle 6 = 90°. Allowing different position angles for each 
Gaussian would severely increase the difficulty of deprojecting 
the 2-D model to a 3-D stellar distribution. We have chosen to 
work with only one position angle. The inclination angle for de¬ 
projection is left as a free parameter within the 88 ° < i < 92° 
interval. 

Both geometries are fitted five times, each time with a dif¬ 
ferent random seed for the genetic algorithms, hence leading to 
a different solution. The variation in the resulting parameters of 
each individual fit is used as an indication of the uncertainty and 
the ability of FitSKIRT to constrain them. The difference be¬ 
tween the models lies in the adopted dust distribution, which we 
discuss below. The results of the fitting are summarised in Table 
3, which lists the parameters of the dust geometries, and Table 
2 , which lists the corresponding dust masses. 


30%. The noise in the residual images is a combination of ob¬ 
servational noise and the Monte Carlo noise from the RT sim¬ 
ulations. This is also why the models (middle column) appear 
noisier than the observations (left column). The dust lane in the 
model does seem to extend up to larger distances, while it falls 
off more abruptly in the observed images. In the central part, 
there is a difference of up to a factor of 2 in luminosity. The 
model slightly overestimates the luminosity in the nucleus of the 
galaxy in favour of a better fit to the outer regions. Since the S/N 
in the centre is highest, these pixels will have the most weight 
in the fitting algorithm. A sufficient number of well-fitting pix¬ 
els with lower S/N are needed to compensate for a deviation in 
the model from the observations in these high-S/N pixels to still 
yield a lower x 2 . NGC 4370 is not perfectly symmetric, as is 
already evident from Fig. 1. 

The western side of the galaxy shows a thicker dust lane with 
stronger extinction than the eastern side. On the other hand, the 
exponential dust disc is an axially symmetric analytical geome¬ 
try. Therefore, it is not surprising that some asymmetry can be 
spotted in the residuals. It appears that the FitSKIRT optimal 
solution favours fitting the western part of the dust lane. As a 
result, the model overestimates the extinction in the eastern part 
by about 50% in the u band and 20% in the z band. The spread 
on the five independent fit results is not more that 10 % except for 
the dust scale height and the dust mass, where it amounts to 50%. 
We find a total dust mass of (5.0±2.4)x 10 6 M o for this model. A 
truncated exponential disc could potentially produce better fits. 
However, it would not change the total dust mass significantly 
because the mass in the outer parts is low due to the exponen¬ 
tial fall off. Furthermore, the central over-density will persist. 
We therefore chose not to pursue a fit with these additional free 
parameters, but explored the possibility of a ring geometry. 

Model B 

Dust could be distributed in a ring-like geometry. This has been 
observed for several ETGs (Goudfrooij et al. 1994a; Bendo et al. 
2006; Patil et al. 2007; Finkelman et al. 2008), but is difficult to 
confirm through observations alone when seen edge-on. We ex¬ 
plored a modified ring distribution, which has a Gaussian profile 
in the radial direction and an exponential fall-off in the vertical 
direction. It is characterised by four free parameters: the radius 
Ro at which the highest density occurs, the radial width W of 
density peak, the vertical scale height 77, and again the dust den¬ 
sity in the centre po. The density in the dust ring is described by 


Model A 


The first geometry consists of a classic exponential dust disc, 
generally applied in the RT models of edge-on spiral galaxies 
(Xilouris et al. 1999; Popescu et al. 2000; Bianchi 2007; Baes 
et al. 2010; Popescu et al. 2011; De Looze et al. 2012a; De 
Geyter et al. 2013, 2014) and also derived from several FIR ob¬ 
servations (see e.g. Munoz-Mateos et al. 2009; Verstappen et al. 
2013). The geometry is determined by three free parameters: the 
radial scale length Hr, the axial scale height h z , and the central 
dust density po. The dust density p^, is then described by 


Pd(R, z) = po exp 


Hr h z 


(9) 


from which the total dust mass can be derived. 

The observed images are reproduced relatively well (see 
Fig. 6 ), with residuals (right column) generally staying below 


p d {R, z) = po exp 


(R-Rq ) 2 

2W 2 


exp 



( 10 ) 


The observed images are also accurately reproduced by the 
dust ring geometry (see Fig. 7). The spread on the five inde¬ 
pendent fit results is slightly greater than for the disc geometry. 
Especially the width W and dust mass have large variations (up 
to a factor of 2.5). Nevertheless, the five independent fits yield 
residual maps that show very similar morphology. There is still 
a slight enhancement visible in the centre of the residuals, which 
is, however, less prominent than in Model A, and the discrepancy 
is not more than 50% in most pixels. This is also an axisymmet- 
ric geometry, and again the extinction in the eastern part of the 
galaxy is overestimated, while the western part is reproduced 
better. Judging from the low residual values in the outer parts of 
the dust lane (30-50%), Model B is more accurate than Model A 
in these regions. This is also visible in the model images, where 
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Fig. 6: FitSKIRT model results for the exponential dust disc geometry. Left column are the observations, with masked-out foreground stars. 
Middle column shows the corresponding model images. The residuals are shown in the right column. 


the dust lane is more compact than Model A and corresponds 
better to the observations. Based on these considerations, a dust 
ring seems more likely than an exponential dust disc in NGC 
4370, although the results lie close to each other. The total dust 
mass for this model is (3.3 ± 2.7) x 10 6 M o which is lower than, 
but consistent with Model A. 


Attenuation maps 

To investigate the validity of the FitSKIRT modelling further, we 
produced the attenuation maps for the ugiz bands. SKIRT offers 
the option of also outputting a transparent model (i.e. only stel¬ 
lar emission), which serves a dust-free model. These were then 
convolved to the same PSF as the observed images and com¬ 
pared to the FitSKIRT best-fitting images using eq. (8). Figure 8 
shows the resulting attenuation maps, together with the ones ob¬ 
tained from the MGE fitting in section 3.3. Attenuation maps 
offset the dusty and dust-free images and consequently enhance 
deviations from a smooth dust distribution. In the case of the 
FitSKIRT attenuation maps, which have smooth distributions as 
an input, this is limited to the Monte Carlo noise inherent to the 


RT simulations. For the MGE modelling, which works with the 
observed images as input, observational noise is added on top of 
the intrinsic deviations of a smooth dust distribution. With this 
in mind, we limit ourselves to a qualitative comparison of the 
attenuation maps. 

The MGE attenuation maps appear thin, but with constant 
thickness. At the edges, the attenuation seems to be slightly 
puffed up in the vertical direction. There is also a clear decrease 
in attenuation going from u to z, which is inherent to the extinc¬ 
tion properties of dust. The FitSKIRT attenuation maps exhibit 
the same decrease with increasing wavelength. There are some 
evident differences between the MGE maps and the FitSKIRT 
maps. In the case of the exponential disc geometry, the attenu¬ 
ation features are more elongated and are thicker in the middle 
than at the edges (as is inherent to the analytical distribution). 
As a result, they overestimate attenuation values in the centre 
and outskirts of the galaxy. The attenuation maps from the ring 
geometry lie closer to the ones from the MGE models. They are 
more compact in scale length and less peaked in the centre. Still, 
the ring geometry yields higher attenuation values in the cen- 
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Fig. 7: FitSKIRT model results for the ring dust geometry. Left column are the observations with masked-out foreground stars. Middle column 
are the corresponding model images. The residuals are shown in the right column. 


tre and fails to reproduce the puffed-up morphology in the outer 
parts. 

In general, the radiative transfer models yield an attenua¬ 
tion map that has higher values (by an average factor of 2), 
but that resembles the observed map. This can partially con¬ 
tribute to the discrepancy in dust mass with the MGE-based es¬ 
timates. However, all models produce roughly the same thin and 
elongated morphology. SKIRT treats the extinction of dust in a 
self-consistent way, so it is reassuring that the attenuation maps 
are reproduced this well. This boosts confidence in the applied 
method and the resulting dust mass parameters. It also under¬ 
lines the importance of the dust geometry and a proper treatment 
of absorption and multiple scattering to convert these attenuation 
values to more reliable dust masses. 

Effect of a stellar disc 

As already mentioned in section 3.2, NGC 4370 may have a nu¬ 
clear stellar disc. This is obscured by the dust lane in the optical 
images, but visible in the IRAC 3.6 /mi image (see Fig. 1). Such 
a component is not explicitly included in the RT models when 


deriving the stellar distribution from the z band. We explored 
the possibility of deriving the stellar distribution from the IRAC 
3.6 jim image using an MGE model, similar to the one derived 
from the z band image. The drawback to this method is the dif¬ 
ference in spatial resolution between IRAC and the NGVS data. 
We only explored the ring geometry and ran the FitSKIRT sim¬ 
ulation five times to get a grip on the uncertainties. 

Figure 9 shows the results of the FitSKIRT simulation. At 
first sight, the observations are reproduced quite well. The dust 
lane is rather compact, and the outer stellar halo is retrieved. 
The inner stellar distribution, however, is much brighter above 
the dust lane in the model. Below the dust lane, it is fainter. The 
larger beam size of the IRAC image also leaves its mark: the 
stellar component is generally more stretched along the major 
and minor axes. The quality of the fit is easier to judge from the 
residuals, as shown in the right-hand column of Fig. 9. There is a 
large residual feature above the dust lane due to the overestima¬ 
tion of the stellar emission. The residuals at the location of the 
dust lane itself again show the features that are present for the 
other models (offsets in the centre and asymmetry in the outer 
part), but they are more prominent here. Just above and below 
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Attenuation 


Fig. 8: Attenuation maps for the various models of NGC 4370. Left: Attenuation maps from the MGE models. Middle: maps from the FitSKIRT 
exponential dust disc model. Right: maps from the FitSKIRT dust ring model. All images are convolved to the same beam size. 


the dust lane, the residuals of the nuclear stellar disc are also vis¬ 
ible. This means the brightness profile of the stellar disc in the 
ugiz bands does not scale linearly with the shape it has at 3.6 jam. 
It is less prominent in the optical bands. This suggests a colour 
gradient in the disc, which is stronger than the general stellar dis¬ 
tribution. Consequently, if the disc is relatively red, the embed¬ 
ded star formation will be low. This would mean that the caveat 
we described at the end of section 3.2 will have a relatively low 
impact on the dust mass estimate from the colour maps or the 
MGE fits. Thanks to the inclusion of a stellar disc, we retrieved 
a higher total dust mass of (7.0± 1. 8 ) x 10 6 M©, which still resem¬ 
bles the other FitSKIRT models. The dust masses of all models 
are consistent with each other within 1 rms. This again suggests 
that the effects of the nuclear stellar disc on the estimate of the 
dust mass are relatively small. 

4. Discussion 

The existence of a diffuse dust component in dust-lane ETGs 
was first proposed by Goudfrooij & de Jong (1995) to explain 
the dust mass discrepancy between optically derived dust masses 


and /RAS-based masses. The foreground dust screen geometry 
they adopted naturally yields lower limits to the optical dust 
mass. In a similar two-component modelling, Kaviraj et al. 
( 2012 ) computed the dust mass in clumpy or filamentary struc¬ 
tures in dusty ETGs. The dust mass in clumps was found to 
be on average five times lower than the /RAS-derived total dust 
masses, again favouring a diffuse dust component, invisible in 
the UV/optical/NIR. A similar discrepancy in optical and IRAS- 
derived dust mass was found for NGC 4370 (Finkelman et al. 
2008). In this work, we re-evaluate this discrepancy using data 
of higher quality at optical and FIR/sub-mm wavelengths. Using 
colour maps, we find an average dust mass of 5. 6 +J 3 x 10 5 M© for 
a foreground screen geometry, which should be seen as a lower 
limit. From MGE modelling, we find Md = 3.4+^ x 10 5 M Q for 
the same geometry, and Md = 7 . 9^5 x 10 5 M© for a uniform 
dust-star mix. On the other hand, modelling of the FIR/sub-mm 
SED yields Mj = (6.2±0.4)x 10 6 M o , and a panchromatic MAG- 
PHYS fit yields Md = 5.7 ^ 6 x 10 6 M©. This quantifies the dust 
mass discrepancy at a factor between 7-18. From the optical 
data, we only found dust inside the dust lane. The spatial resolu- 
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Fig. 9: FitSKIRT model results for the ring dust geometry, but with the stellar geometry derived from the IRAC 3.6 jam image. Left column shows 
the observations with masked-out foreground stars. The middle column gives the corresponding model images. The residuals are shown in the 
right column. 


tion of the FIR/sub-mm data was not high enough to disentangle 
any potential dust at higher scale heights. 

Wise & Silva (1996) argue that a diffusely distributed dust 
component, invisible in single-band optical images, can be 
traced by means of broadband optical colour gradients. Using 
a diffuse dust model, they show that the gradients expected from 
/RAS-based dust masses are compatible with the observed gradi¬ 
ents. More precisely, if the observed colour gradients were com¬ 
pletely due to dust attenuation, the resulting dust masses would 
be significantly higher than the IRAS dust masses. They note 
that age and metallicity gradients may also contribute. There¬ 
fore, their dust masses may be considered as upper limits. In 
the case of NGC 4370, we found that the colour gradients out¬ 
side the dust lane are relatively flat, suggesting only minor dust 
reddening. 

On the other hand, a diffuse dust component is somewhat 
at odds with more recently derived colour gradients (Michard 
2005); models with significant amounts of diffuse dust have 
colour gradients that contradict the observed ones. In fact, such 
objects would completely destroy the correlation they find be¬ 


tween B-R and B-V colours. Additionally, Vanderbeke et al. 
(2011) find a one-to-one correspondence in stellar velocity dis¬ 
persions in the optical and near-infrared. Comparing this to the¬ 
oretical models (Baes & Dejonghe 2000, 2002) suggests there is 
little effect of diffuse dust on the stellar kinematics. They con¬ 
clude that ETGs are virtually optically thin. 

Thanks to the new possibilities offered by FitSKIRT, we are 
able to simulate the dust-starlight interaction properly. The au¬ 
tomated inverse RT fitting makes it possible to determine the 
best parameters in an unbiased way. The only assumption that is 
being made is the choice of dust geometry. We fit two axisym- 
metric dust geometries to mimic the dust distribution in the dust 
lane: an exponential disc and a ring. We do not add a second, 
diffuse component to the model. The exponential disc geome¬ 
try yields a total dust mass of (5.0 ± 2.4) x 10 6 M 0 , and the ring 
geometry yields (3.3 ± 2.7) x 10 6 M©. While the disc produces 
slightly higher dust masses than the ring, both values are close 
and fall within their respective error bars. The uncertainties on 
these masses are in fact the rms of the dust masses of five equiv¬ 
alent simulations, run for each geometry. We found that the x 2 
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Table 3: Overview of the resulting dust geometry parameters and their 
variation from the FitSKIRT modelling of NGC 4370. The obtained 
dust masses are listed in Table 2. 


Parameter 

Unit 

ugiz ± rms 

Exponential disc 

Scale length h R 

pc 

3966 ± 1953 

Scale height h z 

pc 

155 + 17 

Inclination /disc 

deg 

89.27 + 0.10 

Ring (z band) 

Radius Ro 

pc 

1969 + 498 

Width W 

pc 

2130 + 3346 

Height H 

pc 

146 + 14 

Inclination /Ri ng 

deg 

88.83 + 0.59 

Ring (IRAC 3.6 yum) 

Radius Ro 

pc 

1094 + 2085 

Width W 

pc 

4199 + 1195 

Height H 

pc 

299 + 93 

Inclination /Ri ng 

deg 

89.08 + 0.04 


values for each of the simulations lie very close together, but the 
resulting dust masses vary significantly, hence the relatively high 
rms values. Most simulations actually yield a dust mass that is 
higher than the one with the lowest^ 2 , so the uncertainties on 
the dust masses should be seen as possible fluctuations towards 
the high end. With this in mind, we can say that the total dust 
masses from the RT fits are fully consistent with the FIR/sub-mm 
and panchromatic SED fit estimates. 

Additionally, we tested an alternative stellar distribution de¬ 
rived from the IRAC 3.6 yum image to minimize the effects of 
dust obscuration and include the nuclear stellar disc. We only 
explore a ring geometry. The global model quality from this set¬ 
up is lower than the previous set-ups. The main issue seems to 
be the resolution, which causes a slightly different inclination 
angle and a stretched-out stellar profile. The residual features 
at the location of the dust lane are similar to the ones from the 
previous methods, but are generally larger. We find a total dust 
mass of (7.0 ± 1.8) x 10 6 M 0 , which is slightly larger, but con¬ 
sistent with the other RT methods. In this manner, it seems that 
the effect of a nuclear disc on the dust mass determination is mi¬ 
nor. The dust mass from this model is also fully consistent with 
the FIR/sub-mm and panchromatic SED fit estimates. However, 
given the strong residual features, we do not prefer this model to 
the other RT fits. 

We must point out that the dust mass derived from FitSKIRT 
may be overestimated. An exponential disk has its highest den¬ 
sity in the centre, which is unlikely for the case of NGC 4370. 
Second, the density in this geometry only falls to zero at infinite 
radii. Dust extinction is visible in the outskirts of the galaxy in 
the disc model, while this does not seem to be the case in the 
observed images. These two arguments may lead to an overes¬ 
timation of the dust mass. The ring geometry does not suffer 
from the above problems. Inherent to the analytical description 
of the geometry, there is less dust in the centre and almost none 
at large radii. One can therefore assume that the representation 
of the dust geometry by a ring is better. A final argument, ap¬ 
plicable to both geometries, is a line-of-sight effect. NGC 4370 
is an edge-on galaxy, meaning we only see one side of the dust 
lane. It is possible that the dust distribution does not extend to 
the other side of the galaxy. These so-called "arcs” or incomplete 
dust structures have been observed in several dusty ETGs with 
different inclinations (Goudfrooij et al. 1994a; Patil et al. 2007; 


Finkelman et al. 2008). In the extreme case that the other side of 
the galaxy is completely free of dust, this can cause an overesti¬ 
mation of the total dust mass by a factor of 2. This underlines the 
importance of a realistic dust geometry and the proper treatment 
of absorption and multiple, anisotropic scattering to derive dust 
masses from optical data. 

Our models do not assume anything about the emission of the 
dust in NGC 4370. Such energy balance studies have been con¬ 
ducted for several spiral galaxies (e.g. Baes et al. 2010; Popescu 
et al. 2011 ; De Looze et al. 2012a; De Geyter et al. 2015) and for 
the Sombrero galaxy (Ml04, De Looze et al. 2012b). However, 
this requires assumptions on the dust emissivity, sub-resolution 
physics (clumps and filaments) and a thorough knowledge of the 
dust-heating mechanisms: old and young stars, cosmic rays, and 
potentially external X-ray heating (NGC 4370 lies close to mas¬ 
sive ellipticals in the Virgo cluster). All this requires further in¬ 
vestigation and falls beyond the scope of this paper. 

Foreground screen and uniform mix models run short to rep¬ 
resent the 3D intrinsic dust geometry. They should be seen as 
lower limits. Specifically in NGC 4370, RT simulations show 
that the dust is modelled best with an exponential disc or ring 
geometry. The latter produces lower residuals and images that 
are more similar to the observations. The total mass of the dust 
in these geometries is slightly lower, but consistent with Her- 
schel-derived dust masses within 1 rms. Our results suggest that 
it is premature to invoke a diffuse dust component to explain the 
amount of dust that is observed in NGC 4370. 


5. Conclusions 

In this paper, we explored several methods of deriving the to¬ 
tal dust mass in the dust-lane ETG NGC 4370. Table 2 gives 
an overview of the total dust mass for all methods. We briefly 
summarize our findings here. 

- Exploiting the sensitivity and wavelength range of Her- 
schel , we were able to constrain the dust mass very well within 
the limits of the dust model. We find a dust mass that is an order 
of magnitude higher than previous estimates derived from IRAS 
fluxes. 

- We constructed colour maps based on high-quality ugiz 
data and found a lower limit to the dust mass, which lies an order 
of magnitude below the FIR/sub-mm-based estimate. 

- We find a dust mass comparable to the colour-derived mass 
through MGE fits to the stellar distribution of NGC 4370. These 
dust-free MGE models allow the creation of attenuation maps 
when compared to observations. Assuming a simple foreground 
screen geometry for the dust, this yields dust masses comparable 
to the average colour-derived mass. For a uniform mixing of 
dust and stars, we found dust masses that are two to three times 
higher. 

- Since neither of these relatively simple methods yields a 
dust mass that is high enough to match the FIR/sub-mm dust 
mass, we resorted to the more advanced method of radiative 
transfer simulations. Using FitSKIRT, an inverse radiative trans¬ 
fer optimizer around the radiative transfer code SKIRT, we could 
evaluate more realistic dust geometries: an exponential disc and 
a ring. Both geometries resulted in a decent fit, however the 
ring geometry is preferred. The dust masses of both models are 
slightly lower, but fully consistent with the FIR/sub-mm dust 
masses within 1 rms. 

- We checked the influence of a nuclear stellar disc by adopt¬ 
ing the stellar distribution from the IRAC 3.6 jum image instead 
of the z band image. We only explored the ring geometry, but 
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found a best fit model of lower quality than the previous Fit- 
SKIRT models. The dust mass is a fraction higher, but consis¬ 
tent with these models and with the FIR/sub-mm dust masses. 
This suggests that the influence of a nuclear stellar disc is lim¬ 
ited when determining the dust mass in NGC 4370. 

- Within the current models, there is no need to invoke a dif¬ 
fuse dust component for NGC 4370. Instead, we suspect that 
radiative transfer simulations can reproduce the true dust mass 
with only dust in the dust lane, provided an accurate approxima¬ 
tion of the intrinsic dust geometry is made. 

These findings do not imply that the more empirical methods 
(colour maps and MGE fits) in itself are incorrect. If one adopted 
a proper description of a ring geometry and constructed the at¬ 
tenuation - optical depth formula as in equations (3) or (4), these 
methods would yield more reliable dust masses. We therefore 
caution against over-interpreting dust masses and optical depths 
based on optical data alone, when using overly simplistic star¬ 
dust geometries. In fact, for NGC 4370, the difference in total 
dust mass between the simple screen geometry and the 3D radia¬ 
tive transfer fits is a factor of ~10. 

We aim to expand this effort to a set of ten early-type galax¬ 
ies with prominent dust lanes. The objects in this sample have 
all been observed with Herschel as part of the Far-infraRed In¬ 
vestigation of Early-type galaxies with Dust Lanes (FRIEDL). 
This will allow us to check the conclusions of this work against 
similar objects. 
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